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1 .  INTRODUCTION 


In  an  earlier  study  (Rachele  et  al.  ,  1991),  we  concluded  that  to  estimate  turbu¬ 
lence  characteristics  such  as  the  friction  velocity ,  temperature  scaling  length, 
and  the  Obukhov  (1946)  length  from  discretely  measured  data.  It  is  desirable  and 
efficient,  but  not  necessary,  to  have  the  sensors  spaced  geometrically  with 
height.  In  that  study  wc  presented  the  mathematics  necessary  to  show  that  geo¬ 
metric  spacing  Is  not  necessary.  The  mathematical  argument  for  geometric  spacing 
was  not  developed  or  presented.  Furthermore,  we  could  not  find  a  reference  for 
It.  Even  so,  we  considered  the  traditional  geometric  spacing  equations  and  com¬ 
pared  results  from  them  with  our  more  general  formulations  based  on  the  mean 
value  theorem  for  a  few  examples  and  found  that  for  practical  applications  the 
differences  were  not  significant. 

For  most  turbulence  studies  In  the  boundary  layer,  one  desires  single  values  of 
the  turbulence  characteristics  through  a  layer  of  perhaps  20,  30,  or  100  m;  that 
Is,  one  typically  assumes  that  the  characteristic  lengths  are  constant  with 
height.  A  procedure  for  determining  the  total  layer  values  is  to  compute  the 
values  In  the  geometrically  spaced  sublayers  and  then  average  the  sublayer 
values . 

In  this  study  we  consider  the  case  where  the  sensors  are  geometrically  spaced. 
We  then  present  our  mathematical  development  requiring  geometric  spacing  for  a 
single  sublayer  based  on  the  mean  value  theorem  and  assume  that  the  values  for 
the  total  layer  (summation  of  sublayers)  can  be  determined  by  averaging  the 
sublayer  values. 

2.  MATHEMATICAL  DEVELOPMENT 

In  the  development  that  follows,  we  focus  on  the  vertical  gradient  of  wlndspeed 
written  In  similarity  form  for  unstable  and  stable  conditions.  Results  for 
temperature  are  similar.  The  neutral  case  results  from  the  reduced  form  of  the 
unstable  or  stable  case. 

2.1  Unstable  Formulations 

We  write  the  wind  gradient  for  unstable  conditions  as 


dz 


kz  \  ^  LI 


(1) 


where 


V  -  wlndspeed 
U*  -  friction  velocity 
k  -  Karman's  constant 
2  -  height 
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0  -  constant 

L  -  Obukhov  length 
Integrating  equation  (1)  results  in 


(2) 


where 


z„  ~  roughness  length, 

z/L  -  the  Monin  and  Obukhov  (1954)  scaling  ratio 


and 


2  tan-^X  +  — 
2 


(3) 


X  = 


as  shown  by  Paulson  (1970) 


(4) 


Kvaluating  equation  (2)  for  two  geometric  heights  Zj  and  Z2,  where  z^  >  Zj,  we 
oh  tain 


V  -  V  =  ^ 
^  ^  k 


(5) 


Dividing  equation  (5)  by  Zj  -  Zj  gives 


V2  -  v^i  _  cr 


^2  -  -  Zi) 


mil 


-  ♦.  (f )  *  ♦. 


(6) 


The  mean  value  theorem  of  calculus  states  that  there  exists  at  least  one  point 
i ,  where  Zj  <  5  <  Z2  such  that 
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(7) 


where 


Using  equations  (6) ,  (7) ,  and  (8)  we  find  that 


(^2  -  Zl)  1 

1  - 

-  ♦,  1 

'  z  ) 

•f  dr 

(^) 

(8) 


(9) 


Now  suppose  that  a  scaling  constant  €  >  0  exists  such  that 


e 


dV 

dz 


(10) 


where  is  the  geometric  mean  of  Zj  and  Z2,  that  is,  Zg  -  (z^  Z2)^^^. 
Using  equation  (1),  we  can  write  equation  (10)  as 


Solving  equation  (11)  for  6  gives 


e  = 


.  (■  - 

'•  (■  -  nr 


(11) 


(12) 


Substituting  ^  from  equation  (9)  into  equation  (12)  yields 
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or  by  substituting  equation  (14)  into  equation  (16) 


(17) 


To  evaluate  equation  (17),  we  require  an  estimate  of  the  Obukhov  (1946)  length 
L.  We  do  this  by  using  a  basic  definition  of  L,  that  is, 
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(18) 


r  = 

kg  0* 


where  0,  Is  the  potential  temperature  at  a  reference  height  and  where  0*  is  the 
temperature  scaling  length  in  the  potential  temperature  gradient  expression  for 
unstable  conditions,  that  is, 


dz  kz  \ 


1 


(19) 


Using  the  same  mathematical  development  for  0*  as  was  used  for  U*,  we  obtained 


where 


(20) 


(21) 


(22) 


Equations  (17),  (19),  and  (20)  are  solved  iteratively  to  obtain  U*,  0*,  and  L. 

For  neutral  conditions,  that  is,  L  —  <»,  we  also  use  these  three  equations,  but 
equations  (17)  and  (20)  reduce  to 


U* 


(23) 


0.  ^  ^2  (^2  ~  Ql) 

In-^ 


(24) 
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We  conclude  this  section  with  a  brief  discussion  of  the  Richardson  (1920)  number 
that  is  one  measure  of  stability.  By  definition  the  Richardson  number  can  be 


written  as 


R  =  3. 

"  6  {dV/dzf 


Substituting  equations  (1)  and  (19)  into  equation  (25)  gives 


^  (i  - 

e  kz\  ^  Lj 
[(l  - 

(kz)^  [\  ^1 


-1/4  12 


If  7  -  equation  (26)  reduces  to 


However,  from  equation  (18) 


„  _  g  kQ*  z 

'''  ■  e 


_£  .  2 

e,  ■  L 


Tticrerore,  equation  (28)  simplifies  to 


=  i 


On  the  other  hand,  if  we  substitute  equations  (17)  and  (20)  into  equation  (27) 
we  obtain 


(in^  -  (t)  *  ♦»  (^)) 


Equation  (30)  allows  one  to  compute  Rj  for  any  Zi  <  z  <  Z2  given  measured  or 
estimated  values  of  0  and  V. 


Finally,  using  equations  (29)  and  (30)  we  have  an  implicit  expression  for  L,  that 
is, 


(62  -  01) 

-]] 

-  ^1)" 

)  ^ 

(i))l 

(31) 


2.2  Stable  Formulations 

The  stable  case  tracks  the  development  of  the  unstable  case,  except  that  some 
equations  take  a  different  form.  As  an  example  equation  (1)  takes  the  form 
(Dyer,  1974;  Hicks,  1976;  and  Hansen,  1980) 


dV  _  (/• 
dz  kz 


(32) 


Equation  (2)  retains  the  form 


V  = 


(33) 


where 


As  with  the  unstable  case,  when  we  integrate  equation  (33)  for  two  geometric 
heights  7.,  and  7.2,  wiiere  7.2  >  Zj,  we  obtain 


(35) 


Dividing  equation  (35)  by  Z2  -  Z2,  gives  an  equation  of  the  same  form  as  equation 
(6) ,  that  is , 


- 


- 


(36) 


13 


Again,  as  for  the  unstable  case,  we  employ  the  mean  value  theorem  to  obtain 


^2  -  ^  av| 

z,  -  z,  dz 

2  1  r 


(37) 


but ,  where , 


s-  a  •  <i) 


(38) 


From  equations  (36),  (37),  and  (38)  we  obtain 


Assuming  that  an  €  exists  such  that 


(^2  ^l) 

-i) 

\  L 

such  that 

P  ^ 

dz| 

/ 

(39) 


(AO) 


we  obtain,  by  use  of  equation  (32), 


+  5-^) 

,  si) 

L/ 

1 

1 

(Al) 


Solving  equation  (41)  for  €  gives 


1+5-:^ 


e  = 


L  ) 


*  =t) 


(42) 
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Substituting  $  from  equation  (39)  Into  equation  (42)  results  in 


Combining  equations  (37),  (40),  and  (43)  results  in 


dz 


z 


sr 


(^^2  -  V,) 

!i " 
\ 

z  ) 

{lI 

l^)) 

Since 


then 


or  by  substituting  equation  (44)  into  equation  (46) ,  we  obtain 


U* 


which  is  of  the  same  form  as  the  unstable  case, 
from  ♦Jyj  unstable. 


However, 


stable 


(43) 


(44) 


(45) 


(46) 


(47) 


differs 
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To  obtain  the  Obukhov  length  L  for  stable  conditions,  we  again  use  the  basic  form 
L  as 


ere. 


(/*«) 


where 


and 


e*  - 


kfO^  e^) 


/ 


In-^ 

2, 


(^49) 


(50) 


3.  TRADITIONAL  FORMULATIONS  AND  APPROACH 


The  traditional  approach  (see  Hansen,  1980)  starts  with  the  basic  form  of  the 
Richardson  number,  that  is, 


R.  =  M  dQ/dz 
^  ®  (dV/dz)^ 


(51) 


which,  as  shown  in  section  2,  reduces  to 


3  iLil 

0 


(52) 


Using  the  following  approximations  for  0*  and  U* 


0*  = 


U*  = 


(53) 


(54) 


16 


equation  (52)  becomes 


Ri 


3^  (^2  -  Oi) 

0  (V^  -  V^f 


P’or  z  -  Zg  equation  (55)  becomes 

R,  =  -  z.  In  ^ 

'  Oz,  (^2  -  V^f  z^ 


(55) 


(56) 


4 .  COMPARISONS 

In  this  section  we  present  comparisons  of  U* ,  0*,  L,  and  Rj  for  four  sublayers, 
(1/2  -2m),  (1  -  4  m),  (2  -  8  m),  and  (4  -  16  m) ,  using  the  proposed  Rachele/ 
Tunick  (RT)  formulations  and  the  traditional  formulations. 

The  values  of  V  and  0  required  to  evaluate  the  equations  were  obtained  from 
profiles  we  constructed  assuming  that  U*  -  0.4  m/s,  0*  =  -0.5  "C,  L  =  -24.747 

ID,  and  z„  -  0.01  m  for  the  unstable  case,  and  U*  -  0.4  m/s,  0*  -  +0.5  “C,  L  •= 
t24./4  7  m,  and  7.„  ”  O.Ol  m  for  the  stable  case. 

Results  are  presented  in  tables  1  through  8. 


TABLE  1. 

U*.  0*.  L,  and  Ri  FOR  THE  SUBLAYER 
(unstable  case) 

(1/2  -  2  m) 

Rachele /Tunick  (RT) 

Traditional 

U*  (cm/s) 

40.0000 

39.84015 

0*  (“C) 

-.49999 

-.496729 

L  (m) 

-24.74775 

-24.76738 

Ri 

-4.04077  X  10'2 

-4.037568  X  lO'^ 

TABLE  2. 

U*,  0*  ,  L,  and  Ri  FOR  THE  SUBLAYER 
(unstable  case) 

(1  4  m) 

Rachele/Tunick  (RT) 

Traditional 

U*  (cm/s) 

40.0000 

39.879406 

0*  ("O 

-.49999 

-.4984686 

L  (m) 

-24.74775 

-24.67877 

Ri 

-8.08154  X  10‘2 

-8.104131  X  10-2 
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TABLE  3.  U*.  6*.  L.  and  FOR  THE  SUBLAYER  (2  -  8  m) 

(unstable  case) 


Rachele/Tunlck  (RT) 

Traditional 

U*  (cra/s) 

40 . 0000 

39.980267 

e*  CC) 

-0 . 50000 

-.5019645 

L  (m) 

-24.74775 

-24.588487 

Ri 

-0.1616308 

-.1626777 

TABLE  4. 

U*.  6*.  L,  and  R^  FOR  THE  SUBLAYER 
(unstable  case) 

(4  -  16  m) 

Rachele/Tunick  (RT) 

Traditional 

U*  (cm/s) 

39.99999 

40.102282 

0*  (“C) 

-.50000 

-.50595617 

L  (m) 

-24.74775 

-24.509748 

Ri 

-0.3232616 

-.3264007 

TABLE  5. 

U*.  0*.  L,  and  Rj  FOR  THE  SUBLAYER  (1/2  -  2  m) 

(stable  case) 

Rachele/Tunick  (RT) 

Traditional 

U*  (cm/s) 

40.00000 

40.634 

0*  (*C) 

0.50000 

0.50793 

L  (ra) 

24.74775 

25.05282 

Ri 

0 . 04041 

0.03327 

TABLE  6. 

U*.  0*  ,  L,  and  R^  FOR  THE  SUBLAYER 
(stable  case) 

(1  -  4  m) 

Rachele/Tunick  (RT) 

Traditional 

U*  (cm/s) 

40 . 00000 

41.332 

0*  CC) 

0 . 50000 

0.51665 

L  (m) 

24.74775 

25.58276 

Ri 

0.08081 

0.05621 
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TABLE  7.  U*,  6*.  L,  and  FOR  THE  SUBLAYER  (2  -  8  m) 

(stable  case) 


Rachele /Tun Ick  (RTl 

Trad i ( I ona 1 

U*  (cm/s) 

40 . 00000 

42.823 

0*  ("C) 

0 . 50000 

0.53529 

L  (m) 

24.74776 

26.63664 

Ri 

0.16163 

0.08576 

TABLE  8. 

U*.  0*.  L,  and  Rj  FOR  THE  SUBLAYER 
(stable  case) 

(4  -  16  m) 

Rachele /Tunick  (RTl 

Traditional 

U*  (cm/s) 

40.00000 

46.093 

0*  ("C) 

0 . 50000 

0.57616 

L  (m) 

24.74777 

28.87072 

Ri 

0.32326 

0.11616 

S.  SUMMARY  AND  CONCLUSIONS 

This  study  was  performed  to  develop  formulations  that  can  be  used  to  compute 
»'st i mates  of  the  height  derivatives  at  the  geometric  mean  height  of  a  sublayer 
(layer  between  any  two  successive  heights)  and  the  similarity  scaling  constants 
lor  I  lu*  sill) layers  using  the  mean  value  theorem  of  calculus.  As  such,  we  renuire 
that  discrete  vertical  wind  and  temperature  data  be  obtained  from  sensors  that 
arc  geometrically  spaced  with  height. 

A  prior  study  (Rachele  et  al.,  1991)  showed  that  the  height  derivative  at  the 
geometric  mean  height  for  similarity  functional  forms  was  not  generally  equal  to 
the  difference  derivative  of  the  discrete  data. 

Values  of  the  scaling  constants,  Obukhov  length,  and  vertical  profiles  of  wind 
and  temperature  obtained  from  the  formulations  developed  in  this  report  were 
compared  with  values  obtained  from  traditional  equations  showing  minor  differ¬ 
ences  for  the  unstable  case  but  somewhat  larger  differences  for  the  stable  case, 
especially  when  considering  Ri. 
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